Shape Effects of Finite-Size Scaling Functions for Anisotropic 
Three-Dimensional Ising Models 
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The finite-size scaling functions for anisotropic three- 
dimensional Ising models of size L\ x L\ x aL\ (a: anisotropy 
parameter) are studied by Monte Carlo simulations. We 
study the a dependence of finite-size scaling functions of the 
Binder parameter g and the magnetization distribution func- 
tion p[m). We have shown that the finite-size scaling func- 
tions for p(m) at the critical temperature change from a two- 
peak structure to a single-peak one by increasing or decreasing 
a from 1. We also study the finite-size scaling near the critical 
temperature of the layered square-lattice Ising model, when 
the systems have a large two-dimensional anisotropy. We have 
found the three-dimensional and two-dimensional finite-size 
scaling behavior depending on the parameter which is fixed; 
a unified view of 3D and 2D finite-size scaling behavior has 
been obtained for the anisotropic 3D Ising models. 



I. INTRODUCTION 

Finite-size scaling (FSS) has been increasingly impor- 
tant in the study of critical phenomena (j]]^]. This is 
partly due to the progress in the theoretical understand- 
ing of finite-size effects, and partly due to the application 
of FSS in the analysis of simulational results. 

Recently, more attention has been paid to the univer- 
sality of FSS functions [H] for both percolation models 
and Ising models [p|-|fO| . It has been shown that sev- 
eral quantities including distribution functions on various 
lattices have universal FSS functions by choosing appro- 
priate metric factors. It should be noted that universal 
FSS functions depend on boundary conditions and the 
shape of finite systems. 

In the percolation problem, it was considered until a 
quite recent time that there exists only one percolating 
cluster on two-dimensional (2D) lattices at the perco- 
lating threshold. But the importance of the number of 
percolating clusters for anisotropic systems, which was 
pointed out by Hu and Lin B, has captured current in- 
terest |5||il]-|T3|] . Aizenman derived the upper and 
lower bounds of the probability for the appearance of 
n percolating clusters, W n , for 2D percolation at criti- 
cality. His result was later confirmed by Monte Carlo 
simulations Jl2| . By use of conformal field theory, Cardy 



p3[ proposed an exact formula for W n at the percolation 
threshold for the systems with large aspect ratios Jl3| . 
The aspect ratio L1/L2 is an important quantity in the 
FSS functions for anisotropic 2D systems of size L\ x L2. 
The "nonuniversal scaling" of the low-temperature con- 
ductance peak heights for Corbino disks in the quantum 
Hall effect has been discussed in terms of the number of 
the percolating clusters JL4| . 

For the Ising model, the aspect ratio dependence of 
Binder parameter |T5| ] at the critical temperature was 
studied by Kamieniarz and B16tc . Binder and Wang 
discussed anisotropic FSS (T^j. The aspect ratio de- 
pendence of the universal FSS functions for Binder pa- 
rameter and magnetization distribution function was dis- 
cussed by Okabe et al. jlO| for the 2D Ising model with 
tilted boundary conditions. For fixed set of the aspect ra- 
tio and the tilt par ameter, the FSS functions were shown 
to be universal . Quite recently, based on the connec- 
tion between the Ising model and a correlated percola- 
tion model, Tomita et al. Q studied the FSS properties 
of distribution functions for the fraction of lattice sites 
in percolating clusters in subgraphs with n percolating 
clusters, / n (c), and the distribution function for magne- 
tization in subgraphs with n percolating clusters, p n (m). 
They studied the change of the structure of the magneti- 
zation distribution function for the 2D system with large 
aspect ratio in terms of percolating clusters. 

Almost all results above for the shape effects on FSS 
functions are for 2D systems except for some works on 
percolation problem . It is quite important to extend 
these arguments for higher-dimensional systems. Con- 
formal invariance plays a role in 2D systems [ p^[ , but 
is not so powerful for three-dimensional (3D) systems as 
for 2D ones. Moreover, there are two limiting cases for 
anisotropic 3D systems, that is, the one-dimensional (ID) 
limit and the 2D limit. It is interesting to study various 
types of scaling for anisotropic 3D systems. 

In this paper, we study the FSS functions for 
anisotropic 3D systems by Monte Carlo simulations. We 
are concerned with the ferromagnetic Ising model on 
the L\ x L\ x aLi simple cubic lattices with the pe- 
riodic boundary conditions, where a is regarded as the 
anisotropy parameter. Attention is mainly paid to the 
effect of anisotropy, and we study the a dependence of 
the FSS functions for quantities such as the Binder pa- 
rameter and the magnetization distribution function. 
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At the critical temperature for the 3D Ising model, 
we find that FSS functions for the magnetization dis- 
tribution function change from a two-peak structure to 
a single-peak one, when the anisotropy parameter a 
is varied from 1. This behavior is observed both for 
the systems with 2D anisotropy and for those with ID 
anisotropy. When a system has a large 2D anisotropy, 
we find that finite-size systems with fixed anisotropy pa- 
rameter show good FSS behavior as 3D systems near the 
critical temperatures for the layered square-lattice Ising 
models. In contrast, when we fix a number of layers and 
apply the FSS analysis for layered systems, these systems 
are scaled as 2D ones. 

We organize the rest of the paper as follows. In Sec. 2, 
we define the FSS functions and describe the quantities 
we treat in this paper. In Sec. 3, we present our simula- 
tional results. In Sec. 4, our results are summarized and 
discussed. 



simple cubic lattices with different values of L\ and a. 
We use the periodic boundary conditions. 
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II. FINITE-SIZE SCALING 

If a quality Q has a singularity of the form Q(t) ~ 
(t = T — T c ) near the criticality t = 0, then the cor- 
responding quantity Q(L,t) for finite systems with the 
linear size L has a following scaling form, 



Q(i,t)^i-^/(ti 1/,y ), 



(1) 



where v is the correlation-length critical exponent and 
f(x) is the scaling function. Of course, the corrections 
to FSS are not negligible for smaller L. The FSS is also 
applicable to the distribution function of Q. At the crit- 
icality t — 0, the FSS function has a following form, 



P (Q;L,t = 0)~L« J ^F(QL"^) 



(2) 



In this paper we focus on the Binder parameter and 
the magnetization distribution function. The Binder pa- 
rameter |15[ is given by 



(m 4 ) , 
(m 2 ) 2 ' 



(3) 



and serves as a measure of the non-Gaussian nature of 
the distribution function. The FSS functions for these 
quantities are given by 



g(L,t) ~ gitL 1 '"), 
p(m;L,t = 0) ~ L f} ' u p{mL l3 l v ) 1 

where /3 is the magnetization exponent. 



(4) 
(5) 
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FIG. 1. (a) Temperature dependence of g for several lat- 
tices with a=l, 1/2, and 1/4. The temperature is represented 
in unit of J. (b) Plot of g as a function of (T - T c )L 1/v . 

First, we show the results for a < 1, that is, a 2D 
anisotropic case. We plot the temperature dependence of 
the Binder parameter g for a=l, 1/2, and 1/4 of various 
sizes in Fig. 0(a). Error-bars are within the size of marks 
unless specified. We give the FSS plots in Fig. 0(b); g is 
plotted as a function of (T — T c )L 1 / !y , where L is given by 
(Li x L\ x aLi) 1 / 3 . For T c , f3 and v, we use the numeri- 
cally estimated values for the 3D Ising model (^0|, that is, 
T c = 4.5114, /3 = 0.320 and v = 0.625. From now on, we 
represent the temperature in unit of J. From Fig. |l|(b) 
we see that the FSS functions have a dependence in 3D 
systems, which is similar to the case of 2D systems. 



III. RESULTS 



We use the Metropolis Monte Carlo method to sim- 
ulate the ferromagnetic Ising model on L\ x L\ x aL\ 
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FIG. 2. p(m)L- /3/,/ at T = T c as a function of mL w " for 
several lattices with a—1, 1/2, and 1/4. 

We can get more information on the a dependence from 
the magnetization distribution function p(m). In Fig. ||, 
we show the scaling plot of the magnetization distribu- 
tion function p(m) at T = T c for various sizes with differ- 
ent anisotropy parameters (a < 1). From the figure we 
find good FSS behavior and also a large a dependence. 
By decreasing the anisotropy parameter a from 1, the 
FSS functions for the magnetization distribution func- 
tion change from a two-peak structure to a single-peak 
one. 

Next turn to the case of a > 1, that is, a ID anisotropic 
case. We plot the temperature dependence of g for a=l, 
2 and 4 of various sizes in Fig. |3|(a) . The FSS plots are 
given in Fig. ||(b). We see similar behavior of the FSS 
functions for a > 1, that is, g take smaller values if we 
change a from 1. We show the a dependence of the FSS 
functions of the magnetization distribution functions at 
T = T C for a > 1 in Fig. [|. By increasing a from 1, FSS 
functions for p(m) change from a two-peak structure to 
a single-peak one again. 
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FIG. 3. (a) Temperature dependence of g for several lat- 
tices with a—1, 2, and 4. (b) Plot of g as a function of 
(T-TJL 1 ^. 



If we compare the results of Figs. |l|(b) and |3j(b), and 
also those of Figs. || and [|, it seems that there could be a 
set of the anisotropy parameters a < 1 and a > 1 which 
give the same FSS functions. But it needs more careful 
calculations to clarify this universality. 
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FIG. 4. p(m)L- /3/l/ at T = T c as a function of mL w " for 
several lattices with a=l, 2, and 4. 

Let us consider the origin of the a dependence of FSS 
functions. This problem is related to the multiple perco- 
lating clusters. The connection between critical phenom- 
ena of spin models and percolation problems has been 
studied for a long time ]2l]-|23]| . Quite recently, Tomita et 
al. have used the cluster formalism to investigate the 
percolating properties of the 2D Ising model. They elu- 
cidated that the existence of several percolating clusters 
for anisotropic finite systems leads to the change of the 
structure of p(m). That is, the combination of up-spin 
clusters and down-spin clusters gives the contribution of 
m ~ in p(m) for anisotropic case. The a dependence of 
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p(m) for the 3D systems can be understood in the same 
way as the 2D systems. We can apply this argument 
both cases for 2D anisotropy (a < 1) and ID anisotropy 
(a > 1). The value of Binder parameter g at T = T c , 
of course, reflects upon the structure of the distribution 
function p(m). 
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FIG. 5. g at T — T c , g c , as a function of a. We use a 
semi-logarithmic scale. 

In order to see the anisotropy parameter dependence 
more clearly, we plot the a dependence of the Binder 
parameter g at criticality, g c , in Fig. |^. Here, we use 
the semi- logarithmic scale for plot. We add several data 
other than those given in Figs, [l] and |[ For large enough 
a or small enough a, the corrections to FSS become 
larger. From Fig. |^, we find that g c takes the maximum 
value at a = 1 and decreases gradually in both direc- 
tions, a < 1 and a > 1. For very large a or very small 
a, g c tends to vanish, which indicates that the distri- 
bution function p(m) approaches the Gaussian distribu- 
tion. This behavior is consistent with above-mentioned 
multiple-percolating-cluster argument. Such an a depen- 
dence of g at criticality for ID anisotropy of 2D Ising 
model was already discussed for 2D Ising model |Q . 

It is also interesting to consider the relation to the 
layered square-lattice Ising model for the case of large 
2D anisotropy. The critical temperatures and the shift 
exponent for the layered square-lattice Ising model were 
studied by Kitatani et al. p3| . In Fig. ||, we give the 
temperature dependence of g for several lattices, where 
the number of layers, s, is 4, 6, and 8. If we make the 
system size large with fixing s, the system becomes the 
layered square-lattice Ising model. We can estimate the 
s-dependent critical temperature T cs from the crossing 
point of the curves for fixed s in Fig. ^[ Our estimates, 
T c4 = 4.2364(4), T c6 = 4.3701(4) t _and T c8 = 4.4222(4), 
are consistent with those of Ref. 
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FIG. 6. Temperature dependence of g for several systems. 
The numbers of layers are 4, 6, and 8. 
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FIG. 7. Scaling plot of g for fixed a(= 1/16). The scaled 
variable is (T — T CS )L}^ U for each s. The values of T C3 are 
Tca = 4.2364, T c6 = 4.3701, and T c8 = 4.4222. We use the 3D 
exponent, v — 0.625. 

Picking up the data for a = 1/16 from Fig. ^, we plot 
g as a function of (T — T C ^)L X I V in Fig. 0. Here we treat 
T — T cs for each s, and v is chosen as 0.625, that is, the 
3D value. We have good FSS behavior in Fig. \ but this 
FSS plot is different from that given in Fig. [1(b), where 
T — T c has been treated. We should note that T — T cs 
can be rewritten as T — T ca = (T — T c ) + (T c — T cs ), 
and T c — T cs oc s~ A . Here A is the shift exponent and is 
shown to be 1/^3D @]- Thus, for fixed a, we have 3D 
FSS for T — T cs . In this case, the value of g c at T = T cs is 
~ 0.92, which is nothing but the g c value for the isotropic 
2D Ising model. It is different from ~ 0.04 for g c value of 
anisotropic 3D Ising model for a = 1/16, which is given 
in Fig. ||. In Fig- ||, we also show the scaling plot of the 
magnetization distribution function p(m) at T = T cs for 
each s, fixing a as 1/16. We have good FSS behavior 
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using 3D exponents. There are two distinct peaks in 
p(m), which is a typical character of p(m) at T c for the 
isotropic 2D Ising model. In contrast, if we pick up the 
data for fixed s, for example, s = 6, these systems are 
expected to be scaled by using 2D critical exponents. 
Since we fix s(= aL\), the scaled variable is given by 

(T — T CS )L\I V . Using the 2D exponent, z^d = 1, we get 
very good scaling plot, shown in Fig. |[ Of course, the 
critical value g c at T = T c q is ~ 0.92, that is, the g c value 
for the isotropic 2D system. 
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FIG. 8. Scaling plot of p(m) at T = T cs for fixed o(= 1/16); 
p(ra)L~^^ v as a function of mlP'" . The values of T cs are 
T c4 = 4.2364, T c6 = 4.3701, and T c8 = 4.4222. We use the 3D 
exponents, /3 = 0.320 and v = 0.625. 
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IV. SUMMARY AND DISCUSSIONS 

To summarize, we have studied the FSS functions for 
anisotropic 3D Ising model by Monte Carlo simulations. 
The anisotropy parameter (a) dependence of the FSS 
functions for the Binder parameter and the magneti- 
zation distribution function has been investigated. We 
have shown that the magnetization distribution functions 
change from two-peak structures to single-peak ones as 
a increases or decreases from 1. This change of distri- 
bution functions may be attributed to the combination 
of up-spin percolating clusters and down-spin percolating 
clusters for anisotropic systems, which is the same as the 
case of 2D systems p8| . We should also note that the 
distribution function p(m) approaches the Gaussian dis- 
tribution for large anisotropy, which is another indication 
of the multiple-percolation-cluster argument. It is inter- 
esting to calculate the probability for the appearance of 
n percolating clusters, W n , for the anisotropic 3D Ising 
model and to make cluster analysis, which will be left to 
a future study. 

Moreover, we have considered the relation to the lay- 
ered square-lattice Ising models for large 2D anisotropy. 
We have obtained 3D FSS behavior near the critical tem- 
perature of the layered square-lattice Ising models, T cs , 
for fixed a. In contrast, when we fix the number of layers, 
we have shown 2D FSS behavior near T cs . Thus, we have 
obtained a unified view of 3D and 2D FSS behavior for 
the anisotropic 3D Ising models. 

We have studied the shape effects of FSS functions for 
3D Ising models. FSS functions also depend on bound- 
ary conditions. Extension of the 2D results for various 
boundary conditions, for example, tilted boundary con- 
ditions [jl0| , to 3D systems will be interesting. This study 
is now in progress. 
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